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r*j , Abstract 

In this paper, we consider a Gaussian two-way relay channel (GTRC), where two sources exchange messages with each other 
. through a relay. We assume that there is no direct link between sources, and all nodes operate in full-duplex mode. By utilizing 

nested lattice codes for the uplink (i.e., MAC phase), and structured binning for the downlink (i.e., broadcast phase), we propose 
two achievable schemes. Scheme 1 is based on "compute and forward" scheme of [Tj while scheme 2 utilizes two different lattices 
for source nodes based on a three-stage lattice partition chain. We show that scheme 2 can achieve capacity region at the high 
signal-to-noise ratio (SNR). Regardless all channel parameters, the achievable rate of scheme 2 is within 0.2654 bit from the cut-set 
outer bound for user 1. For user 2, the proposed scheme achieves within 0.167 bit from the outer bound if channel coefficient is 
larger than one, and achieves within 0.2658 bit from the outer bound if channel coefficient is smaller than one. Moreover, sum 
rate of the proposed scheme is within 0.334 bits from the sum capacity. These gaps for GTRC are the best gap- to-capacity results 
to date. 
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In this paper, we study a two-way relay channel, where two nodes exchange their messages with each other via a relay. This 
can be considered as, e.g., two mobile users communicate to each other via the access point in a WLAN (2] or a satellite which 
C*~) ' enables data exchange between several earth stations where there is no direct link between the stations. 

Two-way or bi-directional communication between two nodes without a relay was first proposed and studied by Shannon in 
0. In this setup, two nodes want to exchange messages with each other, and act as transmitters and receivers at the same time. 
. For this setup, the capacity region is not known in general and only inner and outer bounds on the capacity region are obtained 
C$ . in the literature. 

The two-way relay channel, also known as the bi-directional relay channel, consists of two nodes communicating with each 
other in a bi-directional manner via a relay. This setup was first introduced in (4]; later studied in |4], 0, ]6] and an approximate 
characterization of the capacity region of the Gaussian case is derived. 

The traditional relaying protocols require four channel uses to exchange the data of two nodes whereas the two-way relaying 
protocol in [4| only needs two phases to achieve bidirectional communication between the two nodes. The first phase is referred 
to as the multiple access (MAC) phase, and the second phase is referred to as the broadcast (BRC) phase. In the MAC phase, 
both nodes transmit their messages to the relay node which decodes them. In the BRC phase, the relay combines the data 
from both nodes and broadcasts the combined data back to both nodes. For this phase, there exist several strategies for the 
processing at the relay node, e.g., an amplify-and-forward (AF) strategy J4), a decode-and-forward (DF) strategy (TJ, or a 
compress-and-forward (CF) strategy ||8]. 

The AF protocol is a simple scheme, which amplifies the signal transmitted from both nodes and retransmit it to them, and 
unlike the DF protocol, no decoding at the relay is performed. In the two-way AF relaying strategy, the signals at the relay are 
actually combined on the symbol level. Due to amplification of noise, its performance degrades at low signal-to-noise ratios 
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(SNR). The two-way DF relaying strategy was proposed in J4], where the relay decodes the received data bits from the both 
nodes. Since the decoded data at the relay can be combined on the symbol level or on the bit level, there has been different data 
combining schemes at the relay for the two-way DF relaying strategy: superposition coding, network coding, and lattice coding 
121 • In the superposition coding scheme, applied in |4), the data from the two nodes are combined on the symbol level, where the 
relay sends the linear sum of the decoded symbols from both nodes. Shortly after the introduction of the two-way relay channel, 
its connection to network coding J9] was observed and investigated. The network coding schemes combine the data from nodes 
on the bit level using the XOR operation, see e.g., ifTOl . IfTTI . Ifl2l . ff3l . lfT4l . lfT31 . Lattice coding uses modulo addition in a 
multi-dimensional space and utilizes nonlinear operations for combining the data. Applying lattice coding in two-way relaying 
systems was considered in, e.g., lfT6l , IfTTI , ffT~8l . In general, as in CF or partial DF relaying strategies, the relay node need not 
to decode the source messages, but only need to pass sufficient information to the destination nodes. 

In this paper, we focus on the Gaussian TRC (GTRC) as shown in Fig. [JJ where nodes 1 and 2 want to communicate with 
each other at rates R\ and i?2, respectively, while a relay node facilitates the communication between them. TRC with a direct 
link between nodes 1 and 2 is studied in (19), |20| . EH . If22ll . l23l . Here, we consider a GTRC without a direct link between 
nodes 1 and 2. Our model is an extension of tfT~8l and it is essentially the same as those considered in 0, l24l . IfTTI . 11251 . |26], 
lf27l . Similar to lfl"8l . IfTTI . we apply the nested lattice coding. For a comprehensive review on lattices and their performance 
analysis, we refer the reader to the study of |f28"1 . |f29l , |f30l . Nested lattice codes have been shown to be capacity achieving for 
the AWGN channel [31), (32), AWGN broadcast channel (32) and the AWGN multiple access channel HI. Song and Devroye 
[33), using list decoding technique, show that lattice codes can achieve the DF rate of single source, single destination Gaussian 
relay channels with one or more relays. It is also shown that the lattice CF scheme, which exploits a lattice Wyner-Ziv binning, 
for the Gaussian relay channel achieves the same rate as the Cover- El Gamal CF achievable rate lf33l . Nazer and Gastpar by 
introducing compute-and-forward scheme obtain significantly higher rates between users than previous results in the literature 
over a relay network fUJ. By utilizing a deterministic approach fl34j and a simple AF or a particular superposition coding strategy 
at the relay, (6) achieves full-duplex GTRC capacity region within 3 bits from the cut-set bound for all values of channel gains. 
A general GTRC is considered in IfTTI . and it is shown that regardless of all channel parameters, the capacity region for each 
user is achievable within | bit and the sum rate is within log | bits from the sum capacity, which was the best gap-to-capacity. 
Communication over a symmetric GTRC when all source and relay nodes have the same transmit powers and noise variances, 
is considered in 11351 . It is shown that lattice codes and lattice decoding at the relay node can achieve capacity region at the 
high SNR. Lim et.al, using noisy network coding for GTRC, show that the achievable rate of their scheme for each user is 
within i bit from the cut-set bound. Il36l . The achievable sum rate in lf36l is within 1 bit from the cut-set bound. 

Here, we apply nested lattice codes at the relay to obtain linear combination of messages. By comparing with the cut-set outer 
bound, we show that the achievable rate region, regardless of all channel parameters, is within 0.2654 bit from the outer bound 
for user 1. For user 2, the proposed scheme achieves within 0.167 bit from the outer bound if the channel coefficient is larger 
than one, and achieves within 0.2658 bit from the outer bound if channel coefficient (g) is smaller than one. The achievable 
sum rate by the proposed scheme is within 0.334 bit from the cut-set bound. Thus, the proposed scheme outperforms IfTTI . 0, 
l36l and the resulting gap is the best gap-to-capacity result to date. 

The remainder of the paper is organized as follows. We present the preliminaries of lattice codes and the channel model in 
Section [TTJ In Section HTlJ we introduce two achievable coding schemes. Scheme 1 is based on compute-and-forward strategy 
while scheme 2 utilizes two different lattices for source nodes based on a three-stage lattice partition chain. In Section II VL we 
analyze the gap between the achievable rate region and the cut-set outer bound for each user. Numerical results and performance 
comparison between two schemes and the outer bound are presented in Section fV] Section fVTJ concludes the paper. 

II. Preliminaries: Lattices and Channel Model 

A. Notations and Channel Model 

Throughout the paper, random variables and their realizations are denoted by capital and small letters, respectively, x stands 
for a vector of length n, (xi,X2, ...,x n ). Also, |.| denotes the Euclidean norm, and all logarithms are with respect to base 2. 
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In this paper, we consider a Gaussian two-way relay channel (GTRC), with two sources that exchange messages through 
a relay. We assume that there is no direct link between the sources and all nodes operate in a full-duplex mode. The system 
model is depicted in Fig. Q] Communication process takes place in two phases: MAC phase BRC phase, which are described 
in the following: 

• MAC phase: In this phase, first the input message to both encoders, Wi, W 2 , are mapped to 

XV = fi(W i ,Yl~ 1 l f°r* = 1,2 

where fi is the encoder function at node i and Y*~ l = ^Y^\ Y^ 2 \ ■ ■ ■ , Y^' j is the set of past channel outputs at 
node i. Wit! 
P > 0, i.e., 



node i. Without loss of generality, we assume that both transmitted sequences are average-power limited to 



1 ™ 

n £—< 

t=i 



X 



w 2 



< P, for i = l,2. (1) 



Both nodes send their signals to the relay. The received signal at the relay at time t is specified by 

where and xlp are the signals transmitted from node 1 and node 2 at time t, respectively, y/g denotes the channel 
gain between node 2 and the relay and all other channel gains are assumed to be one. Z R represents an independent 
identically distributed (i.i.d.) Gaussian random variable with mean zero and variance Nr, which models an additive white 
Gaussian noise (AWGN) at the relay. 

BRC phase: During the broadcast phase, the relay node processes the received signal and retransmits the combined signals 
back to both nodes, i.e., the relay communicates signal Xr to both nodes 1 and 2. Since the relay has no messages of its 
own, the relay signal at time t, X^ , is a function of the past relay inputs, i.e., 

= f* (Yr 1 ) > 

where Jr is a encoding function at the relay and Y^ 1 = ^Y R \ Y R 2 \ • • • , Y R ^ jis a sequence of past relay inputs. The 
power constraint at the relay is given by 



1 

-Ye 

t=i 

The received signal at each node at time t is given by 



W 2 



<Pr, 



Y^ t] = X^ + Zf\ fori = 1,2 

where Xj^ is the transmitted signal from the relay node at time t and Z t - represents an i.i.d AWGN with zero mean and 
variance TVj. 

Node 1, based on the received sequence, Yi, and its own message, W\, makes an estimate of the other message, W2, as 

w 2 = Vi (^1,^1), 

where ipi is a decoding function for at node 1. Decoding at node 2 is performed in a similar way. The average probability of 
error is defined as 

P e = Pr |Wj ^ Wi or W 2 ^W 2 y 

A rate pair (R\, R 2 ) of non-negative real values is achievable if there exists encoding and decoding functions with P e — > as 
n — > 00 1371 . The capacity region of the GTRC is the convex closure of the set of achievable rate pairs R 2 ). 

Note that the model depicted in Fig. [T] is referred to as the symmetric model if all source and relay nodes have the same 
transmit powers, noise variances, i.e., P = Pr and Ni = N 2 = Nr. 
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Figure 1. System Model. Two communication phases for a Gaussian two-way relay channel: MAC phase and Broadcast phase. 



B. Lattice Definitions 

Here, we provide some necessary definitions on lattices and nested lattice codes BTI . 0], ll38l . 

Definition 1. (Lattice): An n-dimensional lattice A is a set of points in Euclidean space R™ such that, if x,y £ A, then 
x + y E A, and if x € A , then —xeA. A lattice A can always be written in terms of a generator matrix G E Z nx ™ as 

A = {x = zG : z G T 1 }, 

where Z represents integers. 

Definition 2. (Quantizer): The nearest neighbor quantizer Q(.) associated with the lattice A is 

Qa(x) = argmin llac — ill . 

Definition 3. (Voronoi Region): The fundamental Voronoi region of a lattice A is set of points in R" closest to the zero 
codeword, i.e., 

V (A) = {xeR n :Q(x)=0}. 
Definition 4. (Moments): a 2 (A) which is called the second moment of lattice A is given by 

1 \\>i a \ \\x\\ 2 dx 

<*M = k v T dx ■ (2) 

" JV(A) 

and the normalized second moment of lattice A is 

G(A) = ^ , = ^1, 

[f V(A) dx]n Vn 

where V = f v ^ dx is the Voronoi region volume, i.e., V = Vol(V). 
Definition 5. (Modulus): The modulo- A operation with respect to lattice A is defined as 

x mod A = x — Q{x), 
that maps x into a point in the fundamental Voronoi region. 
Definition 6. (Quantization Goodness or Rogers-good): A sequence of lattices A( n > C R" is good for mean-squared error 
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(MSE) quantization if 

lim G(A^A = J-. 



2ire 

The sequence is indexed by the lattice dimension n. The existence of such lattices is shown in 

Definition 7. (AWGN channel coding goodness or Poltyrev-good): Let Z be a length-i.i.d Gaussian vector, Z ~ J\f (0, a 2 z I„^. 
The volume-to-noise ratio of a lattice is given by 

,. . (Vol(V)) 2/ " 
a z 

where a\ is chosen such that Pr {Z V} = e and J„ is an n x n identity matrix. A sequence of lattices is A'"' Poltyrev-good 
if 

lim (i ( A (n \ e) = 27re, Ve G (0, 1) 

and, for fixed volume-to-noise ratio greater than 27re, Pi{Z V™} decays exponentially in n . 

Poltyrev showed that sequences of such lattices exist Il40l . The existence of a sequence of lattices A^ n ' which are good in 
both senses (i.e., simultaneously are Poltyrev-good and Rogers-good) has been shown in 

ted 

coarse lattice and as the fine lattice. 

(Nested Lattice Codes): A nested lattice code is the set of all points of a fine lattice A^ that are within the fundamental 



Definition 8. (Nested Lattices): A lattice is said to be nested in lattice A^ if A( n ) C Aj™\ A (n) is referred to as the 
coarse lattice and A^ as the fine lattice. 

(Nested Lattice Codes): A nested lattici 
Voronoi region V of a coarse lattice A^ n ', 



c = {A 1 nv}. 

The rate of a nested lattice code is 

fl = I log|C| = I log ^) 

n n Vol(Vi) 

The existence of nested lattices where the coarse lattice as well as the fine lattice are good in both senses has also been shown 
in HI' ED- An interesting property of these codes is that any integer combinations of transmitted codewords are themselves 
codewords. 

In the following, we present a key property of dithered nested lattice codes. 

Lemma 1. The Crypto Lemma H301 . H311 Let V be a random vector with an arbitrary distribution over W 1 . If D is independent 
of V and uniformly distributed over V, then (V + D) mod A is also independent of V and uniformly distributed over V. 

Proof: See Lemma 2 in H301 . ■ 

III. Nested Lattice Codes 

A. Relay Strategies 

In this Section, we introduce two strategies for processing and transmitting at the relay. In both schemes, we recover a linear 
combination of messages instead of separate recovery of messages at the relay. 

1) Scheme 1: Compute-and-Forward: The compute-and-forward strategy is proposed in [1 j. In this scheme, the goal is recover 
an integer linear combination of codewords, V\ and V2, i.e., we estimate 

Vi+aV 2 , (3) 

where a € Z. Since the transmitted sequences are from lattice codes, it guarantee that any integer linear combination of 
the codewords is a codeword. However, at the receiver, the received signal which is a linear combination of the transmitted 
codewords is no longer integer since the channel coefficients are real (or complex). Also, the received signal is corrupted by 
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noise. As a solution, Nazer and Gastpar 0] propose to scale the received signal by a factor such that the obtained vector is 
made as close as possible to an integer linear combination of the transmitted codewords. 

To reach to this goal, for n large enough, we assume that there exist three lattices , Aj an d A(") such that A(™) C 
A< n) C A 2 n) . A< n) and A 2 n) are both Poltyrev-good and Rogers- good and A^"' is Rogers-good with the second moment 

a 2 (A) = P. 

We denote the Voronoi regions of A[ n \ A [ 2 n) and AW with V x , V 2 and V, respectively. 
Encoding: We choose two codebooks C\ and C 2 , such that 

d = {Ajnv}, 
c 2 = {A 2 nv}. 

Now, for each input node i, the message set {l, 2 nfli } is arbitrarily one-to-one mapped onto Cj. We also define two random 
dither vectors Di ~ Unif (V) for i = 1,2. Dither vectors are independent of each other and also independent of the message 
of each node and the noise. Dither Di is known to both the input nodes and the relay node. To transmit a message, node i 
chooses V 1 E Ci associated with the message and sends 

X t = [V l - Di] mod A. 

Note that by the crypto lemma, X, is uniformly distributed over V and independent of V-j. Thus, the average transmit power 
at node is equal to a 2 (A) = P, so that the power constraint is satisfied. 
Decoding: Upon receiving Y R , the relay node computes 

[aY R + Di + aD 2 ] mod A, 
[aXi + a^/gX2 + ocZr + D 1 + 0D2] mod A, 

[Vx + aV 2 + aXi - (V 1 - D x ) + a^X 2 -a(V 2 - D 2 ) + aZ R ] mod A, 
[Vi + aV 2 + {a - l)Xi + {oty/g - a)X 2 + aZ R ] mod A, 
[T + Z cS ] mod A, 

where (a) follows by adding and subtracting the term V\ +aV 2 , and (b) follows from the distributive law of modulo-A operation 
(29), ID, i.e., 

[[X mod A] + Y] mod A = [X + Y] mod A, 
X,Y E M". The effective noise is given by 

Zcfi — (a - l)Xi + {a^fg - a)X 2 + aZ R , 

and 

T = [Vi + aV 2 ] mod A. 

Since Vi, V 2 , Xi, X 2 are independent, and also independent of Z R , and using Crypto lemma, Z c g is independent of V\ and 
V 2 and thus independent of T. The relay aims to recover T from Y d R instead of recovering Vi and V 2 individually. Due to 
the lattice chain, i.e., A 1 -™' C A^' 1 ' C A 2 n \ T is a point from A 2 . To get an estimate of T, this vector is quantized onto A 2 
modulo the lattice A": 

T = Qa-2 (Y dn ) mod A, 
= Qa 2 (T + Z cff ) mod A, 



(b) 
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where Q\ 2 denotes the nearest neighbor lattice quantizer associated with A2. Thus, the decoding error probability at the relay 
vanishes as n — s- 00 if 

Pr(Z cff i V 2 )^0. (4) 

By using Lemma 8 in |l], we know that the density of Z c g can be upper bounded by the density of Z* ff , which is an i.i.d 
zero-mean Gaussian vector whose variance approaches 



N eq 



((<Vff-a) a + {a-l) 2 )p + a 2 N Rl 
as n — > 00. From Definition [7] this means that Pr (Z* s ^ V2) — > so long as the volume-to-noise ratio satisfies 



( Vol (V 2 )) 2/ ™ 
M(A 2 ) = ^ P^>2ne. 

* * eq 

Therefore, for the volume of each Voronoi region, we have: 

Vol (V,) > {2ireN eq ) n/2 i = 1, 2. (5) 

For the volume of the fundamental Voronoi region of AW, we have: 

p \ n/2 



Vo1{v) = {gW)) ■ (6) 

Now, by using © and (O and definition of the rate of a nested lattice code, we can achieve the following rate for each node: 

^<^ G (A)LaJ ' i = 1 ' 2 " 

Since A is Rogers-good, G (A) 2ire — >• 1 as n — >• 00. Thus, 



< ^ 10 : 2 1 ^ I 



2 



Now, we choose a as the minimum mean-square error (MMSE) coefficient that minimizes the variance of the effective noise 
N eq . Thus, we get 

(a^g + l)P 



MMSE ^ + p + Nr ■ 

By inserting Q in (HJ, we can achieve any rate satisfying 

P 




("MMSEx/ff - fl ) 2 + («MMSE - !) 2 ) P + "umSE^, 



(8) 



(9) 



i.e., it is possible to decode T within arbitrarily low error probability, if the coding rates of the nested lattice codes associated 
with the lattice partition A^™) C A^ and A'") C A 2 satisfy ©. In this scheme, we must only obtain an integer linear 
combination of messages. Since we want to obtain an estimate of V± + aV 2 , a should be an integer. On the other hand, we 
aim to obtain higher achievable rates as much as we can. To reach this goal, we choose a as the closest integer to ^fg i.e., 

a = \y/g\ ■ 

2) Scheme 2: Our proposed scheme: Let us first consider a theorem that is a key to our code construction. 

Theorem 1. /g2]? For any Pi > P 2 > 0, a sequence of n-dimensional lattice partition chains A^/A^/A^, i.e., Aj™ C 
A^ C Aq, exists that satisfies the following properties: 

. A[ n) and A 2 n) are simultaneously Rogers-good and Poltyrev- good while Aq is Poltyrev-good. 
. For any e > , P t - e < a 2 (a^) < P h i £ {1, 2} for sufficiently large n. 
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• The coding rate of the nested lattice code associated with the lattice partition A^/A^ is 

^4 iog ™4 log (Sit)) =7+0n(1) ' 

where C\ = |a^' n Vi| and o n (1) — > as n — > oo. The coding rate of the nested lattice code associated with A^/A^ 
is given by 

*-N*i>-X-E{$)-*+3*(§ 

where C 2 = |a^ i) n V 2 |. 

Proof: The proof of theorem is given in ll42l . ■ 
In the following, by applying a lattice-based coding scheme, we obtain achievable rate region at the relay. Suppose that there 
exist three lattices Aj" , A 2 ™^ and A3 = ^VgA^ 1 , which are Rogers-good (i.e., lim G (A-"M = for i = 1,2,3), and 

n— too V / 

Poltyrev-good with the following second moments 



a 2 (A 1 )=P, and a 2 (A 3 ) = gP, 



and a lattice a[j which is Poltyrev-good with 



A^ C Ag l) , 
A^ C A$\ 



Encoding: To transmit both messages, we construct the following codebooks: 

Ci = {A c nVi}, 
c 2 = |^nv 2 



v/5 

Then node i chooses Vj 6 Cj associated with the message Wj and sends 

Xt = [V., - Di] mod Aj, 

where D\ and D 2 are two independent dithers that are uniformly distributed over Voronoi regions Vi and V 2 , respectively. 
Dithers are known at the source nodes and the relay. Due to the crypto-lemma, Xi is uniformly distributed over V, and 
independent of Vi. Thus, the average transmit power of node i is equal to P, and the power constraint is met. 
Decoding: At the relay node, based on the channel output that is given by 

Y R = X x + ^X 2 + Z R , (10) 



we estimate Vi + ^fgV 2 . Depend on the value of g, we consider two cases: 



Case (I): g > 1: Based on Theorem Q] we can find two lattices, A^ and A3, such that A^ C A^ . With this selection 



of lattices, the relay node performs the following operation: 



Y dR = [aY R + D X + ^/gD 2 ] mod Ai 

= [aXi + a^fgX 2 + aZ R + Di + ^/gD 2 ] mod Ai 

= [V 1 + ^V 2 + aX 1 -(V 1 -D 1 )+a^X 2 -^(V 2 -D 2 ) + aZ R ] mod Ai 

= [[Vi + ^gV 2 ] mod A 3 + (a-l)X 1 + ^g(a-l)X 2 + aZ R ] mod Ai (11) 

= [T + Z cS ] mod Ai, 



January 29, 2013 



DRAFT 



9 



where (c) follows from A3 C A^ and A3 



y/gh^" and the distributive law of modulo-A operation. The effective noise 
{a - 1)X 1 + yfg(a - 1)X 2 + aZ R ] mod A X) 

T = [Vi + ^V 2 ] mod A 3 . 



is given by 



and 



Due to the dithers, the vectors Vi, V 2 , Xi, X 2 are independent, and also independent of Z R . Therefore, Z c g is independent 
of V\ and V 2 . From the crypto-lemma, it follows that T is uniformly distributed over ^/gC 2 and independent of Z e g ll42ll . 

The problem of finding the optimum value for a when the lattice dimension goes to infinity, reduces to obtain the value of 
a that minimizes the effective noise variance. Hence, by minimizing variance of Z e g, we obtain 

{9 + l)P n2 , 

MMSE (g + l)p + N R K ' 

The relay attempts to recover T from Y d R instead of recovering V\ and V 2 individually. The method of decoding is minimum 
Euclidean distance lattice decoding 1311 , (43), l40l , which finds the closest point to Y d R in Aq. Thus, the estimate of T is 
given by, 

T=Q Ac (Y dR ). 

Then, from the type of decoding, the probability of decoding error is given by 

p e = Prjr^r} 

= ?i{Z eS ^V c }. 
Now, we have the following theorem which bounds the error probability. 
Theorem 2. For the described lattice partition chain and any rate Ri satisfying 



1 



P 



,.9+1 N R/ 

the error probability under minimum Euclidean distance lattice decoding is bounded by 

p e = e -"(M 22K ~ Rl) )-°"«), 

where E p (.) is the Poltyrev exponent, which is given by H40V 

§-!(x-l-lnaO, l<x<2 



Ep(x) = 



2 < x < 4 
x > 4 



(13) 



and \x] 



max (0, x). 



Proof: The proof of theorem is similar to the proof of theorem 3 in 114211 and removed here. ■ 
Since Ep (x) > for x > 1, the error probability vanishes as n — > 00 if R\ < R\. Thus, by Theorem Q] and Theorem [2] the 
error probability at the relay node vanishes if 



Ri < 

R-2 



< 



- log 

- log 

2 5 



1 


P 


.9 + 1 




9 




.9 + 1 


N R 



(14) 
(15) 
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Figure 2. Comparison of the achievable rate Ri for different values of g. 



Clearly, using a time sharing argument the following rates can be achieved: 



i?i < u.c.e 



R-2 



< u.c.e 




where u.c.e is the upper convex envelope with respect to . 

At low SNR, i.e., SNR < pure (infinite dimensional) lattice-strategies cannot achieve any positive rates for Ri as shown 
in Fig. [2] Hence, time sharing is required between the point SNR = and SNR*, which is a solution of the following equation: 

«™* 

where f(x) = ^log f^j- +xj- We also evaluate numerically the achievable rates for R\ with lattice strategies for different 
values of g. As we observe, with increasing g, the achievable rate with lattice scheme decreases. As it is shown in Fig. [2] the 
maximum difference between two extreme cases (g = 1 and g = oo) is 0.1218 bit. 

Case (II): g < 1: By using TheoremQ] we can choose two lattices and such that Aj"^ C A3 . The relay calculates 
Y d R = [aY a + Di + y/gZ^] mod A3. The equivalent channel is given by 



Y dn = [aXi + a^/gX-2 + aZ R + Di + y/gD 2 ] mod A 3 

= Wi + ^gV 2 + aX x - (Vi - Di) + a^X 2 -y/g(V 2 - D 2 ) + aZ R ] mod A 3 

= [Vi + ^gV 2 + (a- l)X t + - l)X a + aZ R ] mod A 3 

= [T + Z cff ] mod A 3 , 



(16) 
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where (d) follows from Aj™ C Ag™' and A 3 ™^ = Jgh^' and the distributive law of modulo-A operation. The effective noise 
is given by 

Z off = [{a - l)Xi + ^{a - 1)X 2 + aZ R \ mod A3, 



and 



[Vi + y/gV 2 ] mod Ai. 



Due to the dithers, the vectors Vi, V 2 , Xi, X 2 are independent, and also independent of Zr. Therefore, Z c g is independent 
of V\ and V 2 . From the crypto-lemma, it follows that T is uniformly distributed over C\ and independent of Z c g l42l . In 
order to achieve the maximal rate, the optimal MMSE factor is used, i.e., 

a MMSE (g + l)P + N R ' ( ' 

Similar to case g > 1, instead of recovering V\ and V 2 separately, the relay recovers T. Again, the decoding method is 
minimum Euclidean distance lattice decoding, which finds the closest point to Ya R in A^ 1 '. Thus, the estimate of T is given 

by, 

T = Q\ c (Y dR ). 
Theorem 3. For the described lattice partition chain, if 



1 



log 



9 



gP 



g + 1 n r 



i + 



R 2 < R 2 

the error probability under minimum Euclidean distance lattice decoding is bounded by 

p = e -n(i?p(2 2 («2-«2))_ 0ri(1) ) 

where E p (.) is the Poltyrev exponent given in ( 17 iD . 



Proof: The proof of theorem is similar to the proof of theorem 3 in H42H and removed here. ■ 
Here also the error probability vanishes as n — > 00 if R 2 < R 2 since Ep (x) > for x > 1. Thus, by Theorem [T] and 
Theorem [3] the error probability at the relay vanishes if 



Ri < 



R> 



< 



- log 

- log 

2 & 



1 


p 


.9 + 1 


h Nr~ 


g 




.9 + 1 


Nr 



Clearly, using a time sharing argument the following rates can be achieved: 



i?i < u.c.e 



Ro 



< u.c.e 



log 



- log 

2 & 



1 



.9 + 1 
.9 



P 

Nr~ 



gP_ 



(18) 
(19) 



where u.c.e is the upper convex envelope with respect to Jr. 

For SNR < g (g + ij , pure (infinite dimensional) lattice-strategies cannot achieve any positive rate for R 2 as shown in Fig. [3] 
Hence, time sharing is required between the point SNR = and SNR*, which is a solution of the following: 

where f(x) = | log + gxj . We also evaluate numerically the achievable rate R 2 of lattice strategy for different values of 
g. As we see with decreasing g, the achievable rate with lattice scheme is decreased. 
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Figure 3. Comparison of the achievable symmetric rate for different values of g 



B. Broadcast Phase 

We assume that the relay can recover the linear combination of both messages correctly, i.e., there is no error in the MAC 
phase, T — T. The relay attempts to broadcast a message such that each node can recover the other node's message based on 
both the received signal from the relay node and the available side information at each node, i.e., its own message. For the 
decoding at node 1 and node 2, we can use jointly typical decoding or lattice based scheme. Here, we apply jointly typical 
decoding. We consider scheme 2; decoding for scheme 1 is similar to scheme 2. For scheme 2, we also assume that g > 1. Under 
this assumption, we have R2 > R\. Now, we generate 2 nR2 -sequences with each element i.i.d. according to J\f ~ (0,Pr). 
These sequences form a codebook Cr. We assume that there is a one-to-one correspondence between C2 and Cr. 

Let us denote the relay codeword by Xr (t ^. Based on Y2 = Xr + Z2, node 2 estimates the relay message T as T 2 

if a unique codeword Xr, e £b,2 exists such that (^Xr (t' 2 ^J , are jointly typical, where 

Cr,2 = {Xr (V) : T = [Vi + ^v 2 ] mod A 3 , v 2 € C 2 } . 
Note that \Cr.2 \ = 2 nRl . Now, by using the knowledge of V2 and T 2 , node 2 estimates the message of node 1 as: 

Vi = T 2 - sfgv2 mod A 3 . 
From the argument of random coding and jointly typical decoding ll37l . we get 

Ri < - log [ 1 + — 



(20) 



R2 < I log [ 1 



Pr 



Similarly, at node 1, we get 

Rn < 

2 V iVi . 

Now, we summarize our results for both schemes in the following two theorem: 



(21) 
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Theorem 4. For the Gaussian two-way relay channel, the following rate region is achievable: 

to min I - log I -. r I , - log [ 1 + — I 



I I log 7 I , I log f 1 + § I 



, 2 V(( a MMs E Vff- a ) 2 + K™ SE -i) 2 ) J p + «L ISE ^7 ' 2 V ^ 

where a MMSE = (g+^p+jvrc anc ^ a ^ f ' le closest integer to ^fg i.e., a = \y/g\ ■ 

Proof: The proof of theorem follows from the achievable rate-region of scheme 1 at the relay (0 and achievable rates at 
nodes 1 and 2 (ED- ■ 

Theorem 5. For Gaussian two-way relay channel, the following rate region is achievable: 



Ri < min u.c.e 



R 2 < min u.c.e 



1/1 P 



— log , 

2 & \g+l 



,ilo g (i + §] I, 



2 ° V JV- 



— log 

2 S U + 1 ^fl 



, I log f 1 + ^ ) I . 



2 °\ JVi 



Proof: The proof of theorem follows from the achievable rate-region of scheme 2 at the relay ([T8Tl.([T9li and achievable 
rates at nodes 1 and 2 (|20| |. (l2TT i. ■ 

IV. Outer Bound 

By using the cut-set bound, an outer bound for a TRC can be derived. If a rate pair (i? 1; i? 2 ) is achievable for a general 
TRC, then 

R 1 < mm{I(X 1 ;Y R ,Y 2 \X R ,X 2 ),I(X 1 ,X R ;Y 2 \ > X 2 )}, (26) 
i? 2 < mm{I(X 2 ;Y R ,Y 1 \X R ,X 1 ),I(X 2 ,X R ;Y 1 lX 1 )}, (27) 

where the minimization is over a joint probability p (xi, x 2 ,x R ). If we evaluate the outer bound for the GTRC, all the terms will 
be minimized by the product distribution p (xx, x 2 , x R ) = p (x\) p (x 2 ) p (x R ), where p(x\), p(x 2 ) and p{x R ) are Gaussian 
pdfs with zero means and variances of P\, P 2 and P R , respectively. Therefore, from (l26l and (l27i . one can get 11261 

*i < minQlog(l + £),Il 0g (l + |)), (28) 

to < -inQlo g (l + £),ilo g (l + ^)). (29) 

Although there is a gap between the outer region and the achievable rate-region, we show that the gap vanishes at high SNR 
and hence the capacity region is completely determined in this limit of high SNR. 

Corollary 1. At high SNR, the capacity region of the GTRC, is given by 

to < -inQlog(l + £),Ilog(l + ^))- (l), (30) 

to < --Q io g( i+ C)'5 iog ( i+ t))" o(l) ' (3l) 

where o(l) — > as SNR — > oo. 

Proof: Using ( fTTT ) and ( fTBl with a = 1 and taking Ai and A 3 to be lattices that are Poltyrev-good and Rogers-good with 
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the second moment equals to P and gP, respectively, we get the following achievable rate region at the relay: 

Now, by using d20l l and (f2TT >. we can achieve the following rate region for GTRC: 

Rl < mln (l log (^,I,„ g (l + ^)), ,32) 

R2 < nia Q log (^)4, og ( 1 + ^)). (33) 

By evaluating the outer bound in (|28"T) . and (|29] l for SNR — > oo, we observe that we can approach the outer bound at high SNR 
and the proof is complete. ■ 



and the achievable region in ([24), (123) . It 



A. Calculating the gap 

In the following, we bound the gap between the outer region, given in 
is easy to show that the gap for the symmetric model, i.e., Pr = P and N\ = N2 = Nr = N is always greater than the gap 
in other cases. Thus, we focus on the symmetric model. For the symmetric model, the outer bound on the capacity region is 
given by: 



Ri < I log ( 1 



P 

Nj ' 



i?2 < min ( — log ( 1 



N 



, o log 1 



P 

N 



and the achievable rate region is simplified as: 

Ri < u.c.e 



^log 



1 



P 



i?2 < min u.c.e 



.9 + 1 N 

g , gP 



— log 

2 S 



.9+1 iV 



1 + 



(34) 
(35) 

(36) 
(37) 



To obtain the gap between the outer bound and the achievable region, first we bound the gap for R\ . 
1) Gap for R\: Let us define this gap, denoted by £ (P,g,N), as the following: 

i +) 



£(P,g,N) = hog(l 



P 

N 



log 



1 



.9 + 1 



In the following theorem, we provide an upper bound on £ (P,g,N). 

Theorem 6. Let x*be the solution of the equation x i = log e ( x H — jt ) • For any P,g,N, the gap is bounded by £ (P, g, N) 



? (P, g, N) < max {-log x 



.9 + 1 



5+1 



2In2[ a: * + 5 i I 



ulog 1 



g 

3+1 



(38) 



9+1 



Proof: Let us define a; 



A 



AT - 



£{P,g,N) = -log(l + a:)-u.&e. 



log U 



.9 + 1 



Now, we calculate the upper convex envelope with respect to x. First, we compute the line through the origin which is tangent 
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to ilog (x+^j: 



21n2 



where x*is a solution of the following equation 



log e x 



Therefore, we get 



u.c.e 



■log 



9+1 



1 P 



.9 + 1 



1 N 



i log (x + ^ 

X 

21n2 ( 2; . + _^_) 



X > X* 

0<x <x* 



a) For x > x*: £ (x) is given by 



£ 0) = 2 lo § 



!B + 1 



9+1 



log 1 



9 
9+1 



(39) 



9+1. 



Since £ (x) is decreasing with respect to x, hence £ (x) is maximized for x = x* 
b) For < x < x*: £ (x) is given by 



21n2(x* + ^ I 



The maximum of £ (x) occurs at x m = x* ^j, hence we get 



f(x)<£(x m ) = -log x* + 



9+1 



21n2(y + ^) 



(40) 



Combining d39l and d40b completes the proof. ■ 
2) Gap for R2: To obtain the gap for i? 2 , we consider two cases: 

Case 1: g < 1: In this case, the gap between the outer bound d35l l and the achievable rate in ( f3Tb is defined as the following: 



Vl (P,g,N) = ^l0g[l + ^)- 



N 



Theorem 7. Let x* be the solution of the equation x T 

a -+ ^Xl 



9 , 5^ 



.9+1 W 



log e ( gx H — ) . For any P,g,N, the gap for R2 is bounded by 



m {P,g,N) < max { -logf 



gx 



.9 + 1 



gx 



(9+1) 



Proof: Let x = Then the gap is 



Vi (P, 9, N) = - log (gx + 1) - u.c.e 



21n2( 5 x* + ^ I ) 



^ log [ gx 



» n l0 S 



.9 + 1 



gz* + 1 



(41) 



Here, again we calculate the upper convex envelope. In a similar way, we obtain 



1 



log gx + 



.9+1 



x ■ 

< x < x* 
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where x*is the solution of the following equation 



9+1 



l0g e gx + 



.9+1 



Now, we calculate ?/i (%) as follows: 
a) For x > x*: fji (x) is given by 



m (x) = I log ( ) = \ log ( 1 



l 

9+1 



9+1 



gx 



9+1 . 



Since rji (x) is decreasing with x, hence f)i (x) is maximized at x = x*. 
b) For < x < x*: fji (x) is given by 

1 x 
fji (x) = - log (gx + 1) — 

-= 2^2^* + ^- 

The maximum of rji (x) occurs at x* PfT)' nence ' we § et 



rji (x) < fji I x* - 



9 (.9+1) 



■9+1/ 2hx2[gx* + ^ l 



Case II: g > 1: Let x* be the solution of the equation x ^ ± = log e (gx + ■ Then, if there is no intersection between 
the straight line y = — — j-^ — j-^r and the curve \ log (l + the outer bound in (T34t meet the achievable rate i?2 in 



Thus, the gap is zero. In other wise, we evaluate the gap as follows: 

772 (P,g,N) = ^ 1 °8'( 1 + ^ I 



i n ( 9 qP 

- log — h — 

2 8 \g + 1 N 



(42) 



If SNR is fixed, then 772 (P,g,N) is a decreasing function of g; thus 772 (P,g,N) is maximized at g = 1. In the following 
theorem, we bound 772 (P, g, N) for g = 1. 

Theorem 8. Le? fee f/ie solution of the equation -~r = log e (2; + i). For any P and N, the gap for R2 is bounded by 

t] 2 (P,l,N) < 0.167 

A 



= V2 (x) 



Proof: Let x = ^. We can rewrite (l4Zt as 



»?2 (-P, g,N) = - log (1 + x) - u.c.e 







+ i 


jlog 


(h: 











a) For x > x*: 772 (x) is given by 



Since 772 (a;) is decreasing with respect to x, 772 (x) is maximized at x = x* = 1.655. 
b) For < x < x*: 772 (x) is given by 

^ ( " ) = ^ l0S(1 + X) ' 21n2(^ + l) - 
The maximum of 772 (x) occurs at x* — 4, hence, we get 



m (x) < 772 ( x* - - 



The Theorem holds since 772 (x*) < 772 (x* — 3). 
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3) Maximum gap for R\, R 2 and R\ + R 2 : 

Theorem 9. For g < 1, the outer bound in $34} is within 0.2654 bit from the achievable rate in ( 1561 ) for user 1. In addition, 
the outer bound of ( 1551 ) is within 0.2658 bit from the achievable rate in (1371) for user 2 . 

Proof: The gap £ (P, g, N) in ( f38l > is an increasing function of g, hence this gap is maximized for g = 00. Thus, the 
maximum of gap for Ri is 0.2654 bit. For g < 1, the gap in ( f4TT > is decreasing with respect to g. Therefore, this gap is 
maximized as g — > 0. By evaluating d4TT > for g — > 0, we observe that the maximum gap is 0.2658 bit. ■ 

Theorem 10. For g > 1, the outer bound in ( 1341 ) is within 0.2654 bit from the achievable rate, given (1361) for user 1. For user 
2, the outer bound in ( 135t is within 0.167 bit from the achievable rate in ( 13 7| ). 

Proof: The gap in d38l is an increasing function with respect to g; hence this gap is maximized for g = 00. Thus, the 
maximum of the gap for Ri is 0.2654 bit. Now, by using Theorem|8] we observe that the maximum of the gap for R2 is 0.167 
bit. ■ 

Theorem 11. For g < 1, the sum-rate gap of our proposed scheme (i.e., scheme 2) is within 0.334 of the outer bound. 

Proof: We have 

C(P, <?, N) = (Ri + R 2 )o - (Ri + R%)i = (Rio - R11) + (R20 - R21), 
By using ( f38l > and (HTT) . we get 



GOP, .9, AO < <max<( \ log ( x* 




max 



2 V .9+1 

1, ( *, 9 \ 9V*--m l, ( gy* + i 

-io S [gy + — )- - — f v -log( ) ?h («) 



21n2( 32/ * + ^T)' 2 V^* + ^T 



where x*is the solution of the equation x 1 = log e ( x + — 3-r ) and y* is the solution of the equation ■ v 1 = log e ( gy + 
Now, by evaluating the gap (l43l l. which is an increasing function with respect to g, we observe that the maximum gap is 0.334 
bit (For more detail, see Section IW 

Theorem 12. For g > 1, the sum-rate gap of our proposed scheme (i.e., scheme 2) is within 0.334 of the outer bound. 
Proof: For the sum-rate gap, we have 

( 2 {P,9,N) = (Rt + R 2 ) - (Ri + R 2 )i = {Rio - Ru) + (R20 - R21). 
The first term, using (l38l , ca be calculated as 



Rio - Ru < max { ilog (x* + -^—) — 9+1 ^ log ( 1 + - 9+1 1 ] ) , 

1 2 V .9+1/ 21n2U* + ^ V x 



9+1. 



where x* is a solution of the equation - | x y = log e (^x + j^jj- For evaluating the second term, we consider two cases: 
1) There is no intersection between the straight line y = — — 7-^ — j-^r and the curve h log (l + ^), where x* is the solution 

of the equation x 1 = log e ( gx + -fr ) ■ For this case, the outer bound (l34l meet the achievable rate R 2 in d36i l. Thus, the 
gap for the second term is zero and it is sufficient to calculate the gap for the first term. The gap in (l44t is increasing with g 
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and is maximized for g = oo. With evaluating d44l for g = oo, we observe that the maximum gap is 0.2654 bit. 

2) If there is an intersection between the straight line y = — — j-^ — and the curve | log (l + -£), where x* is the 

solution of the equation - | x t = log e ^gx + j+rji we evaluate the following gap: 



£(P, 3, AT) < max { - log | 



.9+1 



a 

9+1 



2 In 2 




+ £ (P,g,N), 



where 



1, / 5 9^ 
- log — h — 

2 S \g + 1 AT 



The gap is decreasing with g. Thus, it is maximized for g = 1. By evaluating the gap for g = 1, we get 

t(P,g,N)< 0.334. 



V. Discussion and Numerical Result 
First, we assume that ^fg is an integer number. Then, the achievable rate region of scheme 1 can be calculated as: 



Ri < 
i? 2 < 

and the achievable rate region of scheme 2 is given by: 



log 
log 



1 P 

~g~TT + Nb. 
i p 



9 + 1 N R 



Ri 

i?2 



< 



< 



log 
log 



1 



.9+1 
.9 



P 

N 



i + 



-+ g 4 

1 N 



As we see, if g is an integer number and larger than one, the achievable rate region of scheme 2 is larger than that of scheme 
2. But for g = 1, both schemes have the same performance as the proposed scheme in ll35l . For other integer values of g i.e., 
g = 2,3,4, scheme 2 have better performance than scheme 1. This comparison is depicted in FigH]for 9 = 3 and in FigO 
for g = 10. 

In Table U and Table ITU we provide the gap between the achievable rate of user 1 (Ri) and the outer bound, given in d38l l. 
We observe that when channel gain increases, the gap increases and the maximum gap is 0.2654 bit. For g = 10, we compare 
the achievable rate (Ri) of scheme 2 and the outer bound in Fig|6] As SNR increases, the gap decreases. 

In Table fTTTI and for g < 1, we compare the gap between the achievable rate for user 2 (R2) and the outer bound, given in 
(|4i1 i. As we observe, when the channel gain decreases, the gap increases and the maximum gap is 0.2658 bit. 

Table I 

The gap between the outer bound and the achievable rate FOR USER 1, i?i FOR g > 1 



Channel Gain (g) 


1 


4 


9 


16 


25 


64 


100 


inf 


Gap (bits) 


0.167 


0.2361 


0.2497 


0.2547 


0.2573 


0.2621 


0.2637 


0.2654 
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Figure 4. Comparison of the achievable rates of scheme 1 and scheme 2 for user 1, R2, in the symmetric GTRC. The channel gain is g = 3. 
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Figure 5. Comparison of the achievable rates of scheme 1 and scheme 2 for user 2, R2, in the symmetric GTRC. The channel gain is g = 10. 
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Figure 6. Comparison of the achievable rate of scheme 2 for R± and the outer bound in the symmetric GTRC. The channel gain is g = 10. 



Table II 

The gap between the outer bound and the achievable rate FOR USER 1, i?i FOR g < 1 



Channel Gain (g) 


1 


0.7 


0.6 


0.5 


0.3 


0.1 


0.001 


0.0001 


Gap (bits) 


0.167 


0.146 


0.136 


0.1252 


0.09497 


0.045 


0.000686 


0.0001611 



Table III 

THE GAP BETWEEN THE OUTER BOUND AND THE ACHIEVABLE RATE FOR USER 2, R 2 



Channel Gain (g) 


1 


0.7 


0.6 


0.5 


0.3 


0.1 


0.001 


0.0001 


Gap (bits) 


0.167 


0.1872 


0.195 


0.2038 


0.224 


0.2498 


0.2651 


0.2658 



VI. Conclusion 

In this paper, a Gaussian two way relay channel (GTRC) is considered. By using nested lattice-based coding scheme, we 
obtain two achievable rate regions for this channel and characterize the theoretical gap between the achievable rate region and 
the outer bound. We show that the gap is less than 0.2658 bit for each user, which is best gap-to-capacity result to date. We 
also show that, at high SNR, the lattice based coding scheme can achieve the capacity region. 
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